January 1998 



Detection of Center Vortices in the 
Lattice Yang-Mills Vacuum 



L. Del Debbio a , M. Fabe? , J. Giedt c , J. Greensite d , and S. Olejmk 6 



a Dept. of Physics and Astronomy, Univ. of Southampton, 
Southampton S017 IB J, UK. E-mail: ldd@hepl.phys.soton.ac.uk 

b Inst, fur Kernphysik, Technische Universitat Wien, 
A-1040 Vienna, Austria. E-mail: f aber@kph . tuwien . ac . at 

c Physics Department, University of California, Berkeley, 
CA 94720 USA. E-mail: giedt@socrates.berkeley.edu 

d The Niels Bohr Institute, DK-2100 Copenhagen 0, 
Denmark. E-mail: greensite@nbivms.nbi.dk 

e Institute of Physics, Slovak Academy of Sciences, 
SK-842 28 Bratislava, Slovakia. E-mail: fyziolej@savba.sk 



Abstract 



We discuss the implementation of the "direct" maximal center gauge (a gauge which 
maximizes the lattice average of the squared-modulus of the trace of link variables), and 
its use in identifying Z 2 center vortices in Yang-Mills vacuum configurations generated by 
lattice Monte Carlo. We find that center vortices identified in the vacuum state account for 
the full asymptotic string tension. Scaling of vortex densities with lattice coupling, change 
in vortex size with cooling, and sensitivity to Gribov copies is discussed. Preliminary 
evidence is presented, on small lattices, for center dominance in SU (3) lattice gauge theory. 



1 



1 Introduction 



In a number of recent articles and conference proceedings |]-[3[ we have presented numerical 
evidence in favor of the Center Vortex theory of confinement, which was put forward in 
the late 1970's (§-§[. Our most important tool is the use of the maximal center gauge, 
which, combined with "center projection," allows us to identify the locations of center 
vortices in thermalized lattice gauge-field configurations. It is found that these vortices, 
by themselves, account for the entire asymptotic string tension ("center dominance"). We 
have also found evidence that the monopoles of the maximum abelian gauge lie along 
the center vortices in a monopole-antimonopole chain, and that their non-abelian field 
strength, above the vacuum average, is almost entirely oriented in the vortex direction. 
This opens the way to explain abelian dominance, and monopole condensation, in terms 
of more fundamental underlying vortex configurations. Finally, in ref . || , we have argued 
that the Casimir scaling of higher representation string tensions, formerly a very strong 
argument against the center vortex theory, can in fact be understood in terms of center 
vortices. Independent arguments in favor of the center vortex theory have been presented 



by Kovacs and Tomboulis [10|, who follow a rather different approach but reach similar 
conclusions. 

In this article we will explain, in section 2, the actual implementation of the "direct" 
maximal center gauge, which underlies much of our work, and review the evidence, in the 
direct gauge, that center vortices are responsible for quark confinement.^ We go on to 
study (section 3) how vortex configurations are affected by cooling the lattice, and then 
take up issues related to Gribov copies (section 4). Some preliminary evidence for center 
dominance in 577(3) lattice gauge theory is presented in section 5, followed by a summary 
of our results in section 6. 



2 Center Vortices and Confinement 

Our procedure for locating center vortices in thermalized lattice configurations was inspired 
by earlier work on abelian projection in maximal abelian gauge fTT|| . The idea is to fix a 
gauge (the maximal center gauge) which, in the case of SU (2) gauge theory, reduces the 
full SU(2) gauge symmetry to the center subgroup Z 2 . "Center projection" is a mapping 
of the full gauge field configuration U^(x) onto a configuration Z^x) = ±1, transforming 
as a Z 2 gauge field under the residual Z 2 symmetry. The excitations of a Z 2 gauge field 
are (thin) center vortices, and these are used to locate thick center vortices in the full, 
unprojected gauge- field configuration , clS explained below. 

In fact we have introduced two versions of the maximal center gauge; an "indirect" 
version, in refs. |1, 12 1, and a "direct" maximal center gauge in ref. ||. The indirect gauge 



is a further gauge-fixing within maximal abelian gauge, which reduces the residual U(l) 
gauge symmetry to Z 2 . One begins by fixing to maximal abelian gauge, defined as the 

^ome of this evidence was obtained previously in a slightly different, "indirect", version of maximal 
center gauge 
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gauge which maximizes 

Q = J2^K(x)a 3 Ul(x)a 3 ] (1) 

X,fl 

leaving a residual U(l) symmetry. This gauge makes the link variables as diagonal as 
possible. Abelian link variables are extracted from the diagonal elements of the full link 
variables, rescaled to restore unitarity: 

_ diag[£/n, U22] 



y/UuUh + U 22 U* 22 

(2) 



e 



These variables transform under the residual symmetry like U(l) gauge field link variables. 
The indirect maximal center gauge uses this residual symmetry to bring the A-link variables 
as close as possible to SU(2) center elements ±1, by maximizing the lattice average of 
cos 2 (#). Center projection is then achieved by identifying 

Zfj,(x) = sign [cos (6^ (x))] (3) 

The gauge is "indirect" in the sense that the center is maximized in the abelian link variable 
A, rather than directly in the full link variable U. String tensions can be extracted from 
the center-projected configurations and, although agreement with the asymptotic string 
tension of the full configurations is not too bad, significantly better results are obtained in 
the direct maximal center gauge, as will be seen below. 

The direct maximal center gauge, in SU(2) gauge theory, is defined as the gauge which 
brings the full link variables U as close as possible to the center elements ±7, by maximizing 
the quantity 

R=j:nu,(x)} 2 (4) 

with center projection defined by 

Z^x) = sign[Tr^(x)] (5) 

Again, the projected Zj^x) field transforms as a gauge field under the residual Z 2 symmetry. 
Before going on to discuss numerical results obtained in this gauge, we must first discuss 
how to implement it. 



2.1 Fixing to direct maximal center gauge 

The gauge- fixing is accomplished by over-relaxation [Tj|. Beginning with a thermalized 
but non-gauge-fixed lattice, we sweep through the lattice site by site. At each site x, one 
needs to find the gauge transformation g which maximizes the local quantity 

R * = l {E^^Wl' + E^^^-A)^^)] 2 } (6) 
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Denote 

g(x) = gJ-ig-B 

Ufj,{x) = rf 4 (/i)J + id(fi) ■ a 

U^x-jl) = d 4 (n + A)I - id(fi + 4) • a (7) 

where the choice of signs in front of the various terms proportional to % is made for conve- 
nience. Then 

! 8 / 4 \ 2 

^ = oE E^4(/) (8) 
z i=i \k=i ) 

We have to maximize this quantity subject to the constraint that g is unitary. To this end, 
introduce a Lagrange multiplier 

R = R x + \\(i-jr^ (9) 

Then the conditions for a maximum satisfying the constraint, obtained by differentiating 
R, are 

4 8 
3=1 1=1 

E»2 = i (io) 
k=i 

This can be written as an eigenvalue equation 

DG = \G (11) 

where 

A, = E^K(0 (12) 
i=i 

and the unitarity constraint is the normalization condition 

G-G = l (13) 

At this point, the problem of finding g boils down to finding the eigenvectors of a 4 x 4 
real symmetric matrix, which is achieved by standard methods. There are four eigenvec- 
tors corresponding to four stationary points. The eigenvector with the largest eigenvalue 
corresponds to the gauge transformation g maximizing R x at site x. 

The next step is to apply the over-relaxation algorithm. We transform the links touching 
site x not by g(x) but by g u {x), where 

g = gj-ig-a 

= cos(0)7 — in ■ <rsin(0) 
g u = cos(uj(f))I — in ■ a sm(uj(f)) (14) 
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and we use lo = 1.7 [fnj. This procedure is applied at each site of the lattice, sweeping 
through the lattice several hundred times. The algorithm stops when a convergence crite- 
rion is satisfied. Our criterion is that the lattice average of (|Trf/) 2 changes by less than 
0.00015 after 50 gauge-fixing sweeps. 

Finally, it must be noted that maximal center gauge, like the Coulomb, Landau, and 
maximal abelian gauges, is afflicted with Gribov copies. To alleviate the problem, we make 
several gauge copies of the original configuration by applying random gauge transforma- 
tions, and then gauge fix each copy. The gauge-fixed copy with the largest value of R is 
chosen for data-taking. In practice, three gauge copies seems sufficient, in the sense that 
we don't improve the average value of R very much (or change the final results) by making 
more copies. 

All of the data reported below was obtained in the direct maximal center gauge. A 
portion of those results, displayed in Figures |l|-|^ below, are similar to results obtained 
previously in the indirect version of maximal center gauge [|TJ. 

2.2 Projection vortices locate center vortices 

Having fixed to the direct maximal center gauge by this procedure, we obtain the corre- 
sponding center-projected configuration (a Z 2 lattice gauge field), from (|). The excitations 
of any Z 2 gauge field are line-like (D=3) or surface-like (D=4) vortices on the dual lat- 
tice. We refer to these excitations, in the center-projected configurations, as "projection 
vortices" or just "P-vortices." The question is whether P-vortices in the projected config- 
uration are in any way related to the existence of center vortices in the full, unprojected 
lattice configuration. 

In order to study this question, we introduce the concept of vortex-limited Wilson loops 
W n (x). We say that a plaquette is pierced by a P- vortex if, upon going to maximal center 
gauge and center-projecting, the projected plaquette has the value —1. Likewise, a given 
lattice surface is pierced by n P-vortices if n plaquettes of the surface are pierced by P- 
vortices. As a Monte Carlo simulation proceeds, the number of P-vortices piercing any 
given surface will vary. Define W n (C) to be the Wilson loop evaluated on a sub-ensemble 
of configurations, selected such that precisely n P-vortices, in the corresponding center- 
projected configurations, pierce the minimal area of the loop. It should be emphasized here 
that the center projection is used only to select the data set. The Wilson loops themselves 
are evaluated using the full, unprojected link variables. In practice, to compute W n {C), 
the procedure is to generate thermalized lattice configurations by the usual Monte Carlo 
algorithm, and fix to maximal center gauge as described above. For each independent 
configuration one then examines each rectangular loop on the lattice of a given size; those 
with n P-vortices piercing the loop are evaluated, the others are skipped. 

The test for whether (thin) P-vortices in the projected configuration correspond to 
(thick) center vortices in the full, unprojected SU (2) gauge-field configuration is whether 
the behavior 



5 



is found in the limit of large loop area. The reasoning behind this test has been given 
elsewhere fl],@], but for completeness we repeat the argument here. 

Vortices are created by discontinuous gauge transformations. Suppose loop C, para- 
metrized by x M (r), r G [0, 1], encircles n vortices. At the point of discontinuity 

g(x(0)) = (-l) n g(x(l)) (16) 

The corresponding vector potential, in the neighborhood of loop C can be decomposed as 

A<£\x) = g- 1 5A^(x)g + ig~%g (17) 

with the inhomogenous term dropped at the point of discontinuity. Then 

W n {C) = < Trexp[i jdx^A^] > 

= (-1) 71 < Tr exp[i i dx^SA^} > (18) 

In the region of the loop C, the vortex background looks locally like a gauge transformation. 
If all other fluctuations SA^f 1 are basically short-range, then they should be oblivious, in 
the neighborhood of the loop C, to the presence or absence of vortices in the middle of the 
loop. In that case, if we have correctly identified the vortex contribution, then 

< Trexp[z I da^SA^] > « < Trexp[z j> dx»5Af] > (19) 

for sufficiently large loops, and eq. ([15]) follows immediately. The validity of eq. fllBT) then 
constitutes a test of whether P-vortices, which we are using to select the subscript n of 
W n (C), actually locate center vortices in the unprojected configurations. 

Figure [I] shows the ratio Wi/W , with the single P-vortex associated with W\ located 
at (or touching) the center of the loop. Likewise, Fig. [| shows the ratio W2/W0, with the 
two P-vortices for W% located near the center of the loop. Both figures were obtained from 
a simulation on a 14 4 lattice at j3 = 2.3 (400 configurations separated by 100 sweeps), and 
both appear to be quite consistent with the limiting behavior (|i~5|) .pl 

For very large loops, the fraction of configurations in which no vortex pierces the loop 
(the subensemble used to compute Wq) becomes very small. So as a further check, using 
all the configurations, we define W evn (C) to be the Wilson loop evaluated in sub-ensemble 
in which only even (including zero) numbers of P-vortices pierce the minimal area, while 
W dd(C) is the corresponding quantity for odd numbers of P-vortices. For a very large loop, 
the fraction of configurations used to evaluate W evn (C), denoted P evn (C), and the fraction 
PoddiC) used to evaluate W a dd(C), should each approach 50% of the total configurations. 
This is in fact the case, as seen in Fig. [| If P-vortices in the projected lattice are associated 
with center vortices in the unprojected lattice, then we would expect, by the same argument 
leading to eq. (|T5D, that 

2 Qualitatively similar results were found in the indirect maximal center gauge H. 
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Figure 1: Ratio of the 1- Vortex to the 0- Vortex Wilson loops, Wi(C)/W (C), vs. loop area 
at = 2.3. 



W oM (C) -> -W em {C) (20) 

in the limit of large loop area. That also appears to be the case, as seen in Fig. |j. 

The conclusion is that P-vortices in center-projected lattice configurations obtained in 
direct maximal center gauge serve to locate thick center vortices in the full, unprojected, 
lattice gauge field configuration. It is well to bear in mind, however, that we have no real 
understanding of why this technique finds center vortices; our confidence is based entirely 
on the numerical results shown in this and the following sections. 

2.3 No vortices means no area law 

The fact that center vortices can be identified in the gauge-field vacuum does not necessarily 
imply that vortices are important for the confinement mechanism. There is, however, a 
simple test of their relevance. Let us define Xn(I, J) as the Creutz ratio extracted from 
the vortex-limited loops W n (C). If the presence or absence of center vortices crossing the 
minimal spanning surface of a loop is unrelated to the area-law falloff, then we would 
naturally expect, at least for large loops, that 

Xo(I,J)*x(I,J) (21) 

where ^(/, J) is the usual Creutz ratio with no restriction on numbers of vortices. In fact, 
the above equation is entirely wrong, as seen in Fig. [5| When Wilson loops are evaluated 
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Figure 2: Ratio of the 2- Vortex to the 0- Vortex Wilson loops, W 2 (C)/W (C), vs. loop area 
at = 2.3. 




Figure 3: Fraction of link configurations containing even (or zero)/odd numbers of P- 
vortices, at (3 — 2.3, piercing loops of various areas. 
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Even, Odd, and Full Wilson Loops 
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Figure 4: Wilson loops W evn (C), W ^(C), and W{C) at larger loop areas, extracted from 
configurations with even (or zero) numbers of P-vortices, odd numbers of P-vortices, any 
any number of P-vortices, respectively, piercing the minimal loop area, again at (3 = 2.3. 



in subensembles in which no vortices cross the minimal area of the loop, the string tension 
vanishes. 

As a further test, we may consider loops pierced by even (or zero) numbers of vortices. 
According to the center vortex theory, the asymptotic string tension is entirely due to fluc- 
tuations in the number of center vortices piercing the surface of the loop. The asymptotic 
effect of creating a vortex piercing (once or an odd number of times) the loop surface is to 
multiply the loop by a center element, i.e. in SU(2) 

Tr[UU...U] -> (-1) x Tt[UU...U] (22) 

and therefore the area-law falloff is due, asymptotically, to a delicate cancellation between 
configurations with even (or zero) numbers of vortices piercing the loop (which gives a 
positive average contribution), and configurations with odd numbers of vortices piercing 
the loop (which gives a negative average contribution). In fact, we have already seen some 
evidence of this cancellation in Fig. f|, where the value of the full loop W(C) is much 
smaller than the magnitudes of either the even or the odd components W evn (C), W dd{C). 
More quantitatively, if we evaluate Creutz ratios x evn {l, J) evaluated from the even- vortex 
(W evn ) contribution alone, the vortex theory predicts that 

Xevn(I,J)^0 (23) 
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Figure 5: Creutz ratios Xo(R, R) extracted from loops with no P- vortices, as compared to 
the usual Creutz ratios x{R-> R) &t (3 = 2.3. 

in the limit of large loop area. Once again, from Fig. |6], this appears to be exactly what 
happens. 

2.4 Center vortices generate the full asymptotic 
string tension 

From the previous results, we deduce that the confining properties of lattice gauge-field 
configurations are strongly correlated with distribution of center vortices. The final check 
is whether these vortices account for the entire asymptotic string tension, as predicted by 
the center vortex theory. We have already seen that the asymptotic effect of n vortices 
piercing the middle of a large loop is to contribute a factor (— l) n to the loop value. In that 
case, the expectation value of a large Wilson loop can be factored into two components: 
(i) a factor W vor (C) due to the effect of vortices crossing the minimal area, far from the 
perimeter of the loop; and (ii) a factor W per (C) due to short-range fluctuations (denoted 
5A^ in eq. fll7D ) around the vortex background, near the loop perimeter. Asymptotically, 
for large-area loops, the vortex theory predicts that 

W(C) -> W vor (C)W per (C) 
W vor (C) = < (-1) 71 > n = no. of vortices piercing C (24) 

Since W per (C) should behave asymptotically like W evn {C) or |Wodd(C)|, it does not have 
an area-law falloff, and the entire string tension must be due to W vor (C). But if, as we 
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Figure 6: Creutz ratios Xevn(R, R) extracted from loops W evn (C), taken from configurations 
with even (or zero) numbers of P-vortices piercing the loop. The standard Creutz ratios 
at this coupling ((3 = 2.3) are also shown. 

have seen in the previous section, P-vortices locate center vortices, then 

W vor (C) =< ZZZ...Z > (center-projected loop value) (25) 

where the product of links ZZ...Z on the projected lattice is taken around loop C. There- 
fore, 

• if P-vortices locate center vortices, and 

• if the center vortex theory is correct; 

• then the string tension of center projected loops should exactly match the asymptotic 
string tension of the full theory. 

Figure |7] is a Creutz ratio plot, extracted from center-projected Wilson loops (i.e. from 
loops of the Zp[x) link variables) in direct maximal center gauge. The straight solid line 
is the usual two-loop expression (cr = string tension, a = lattice spacing) 



2 a /6vr 2 \ 102/121 

aa= TATi p exp 



6n 



2 
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(26) 



with cr/v^A = 58. There are two aspects of this plot which are worth noting in particular. 
First, unlike a standard plot in the unprojected theory, the Creutz ratios almost fall on 
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top of one another, starting at R = 2. This is not so surprising, from the point of view of 
the vortex theory. The short range gluonic fluctuations which give rise to the Coulombic 
potential have been eliminated (these would contribute to W per (C)); only the fluctuations 
in vortex number, which give rise to a linear potential, remain. Second, even 1), which 
is just the logarithm of the center-projected plaquette, appears to be scaling. This fact, as 
we will see, is related to the scaling of the vortex density. 
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Figure 7: Creutz ratios from center-projected lattice configurations, in the direct maximal 
center gauge. 

Scaling of the center-projected string tension is not sufficient for our purposes; what 
is necessary is that the actual value of the string tension, at every /3, agrees with the 
value for the asymptotic string tension of the unprojected configurations. This is also 
what we find. Figure [8] shows our data (triangles) in the scaling region, for Creutz ratios 
x(R, R) of center-projected Wilson loops, as compared to the value for the full theory of 
the asymptotic string tension.0 The values for the full theory (solid lines), with associated 
error bars (dashed lines), are taken from Bali et al. |1^]. This agreement of the center- 

3 This data was taken, for (3 = 2.3 and = 2.4, on 16 4 lattices with 30 configurations separated by 100 
sweeps. A 22 4 lattice and 20 configurations separated by 100 sweeps was used for (3 = 2.5. 
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projected and full asymptotic string tension persists into the strong coupling regime (see 
Fig. below). 



Direct Maximal Center Gauge, /S =2.3, 2.4, 2. 








0.01 



(S = 2.3 



= 2.4 



Figure 8: Center-projection Creutz ratios x(Ri R) vs - R at /3 = 2.3, 2.4, 2.5. Triangles 
are our data points. The solid line shows the value (at each (3) of the asymptotic string 
tension of the unprojected configurations, and the dashed lines the associated error bars, 
quoted in ref. Ui~4"|. 



2.5 Scaling of the center vortex density 

Finally, we consider the density of vortices. In the vortex theory of confinement, vortices 
must be "condensed" in the sense that the average extension of a vortex is on the order 
of the lattice size itself. Its easy to see why. Suppose the opposite were true, i.e. that 
there were some upper limit to vortex extension, and that almost all vortices, in a very 
large lattice, would fit inside a hypercube of side length L. Now consider R x T Wilson 
loops with R,T >> L. Then only vortices within a distance L of the loop perimeter could 
be linked to the loop, and this would lead asymptotically to a perimeter-law, rather than 
area-law, falloff. 



13 



If vortices are physical objects, their density should scale with (5 in some appropriate 
way. P-vortices are located somewhere near the middle of "thick" center vortices in the 
unprojected lattice, and these P-vortices have the topology of surfaces in D=4 dimensions. 
If center vortices scale correctly, P-vortices should also scale. The proper asymptotic 
scaling of P-vortex densities in the indirect maximal center gauge was reported recently by 



Langfeld et al. [0. We can also observe this scaling in the direct center gauge in a rather 
simple way: We first define p to be the fraction, and N vor to be the total number, of center 
projected plaquettes with value —1. N vor is also the total area of all P-vortices on the dual 
lattice, and we denote by Nt the total number of all plaquettes on the lattice. Then 



V 



N N a 2 

ly vor ly vor Uj 



-a 



N T N T a 4 
Total Vortex Area 
6 x Total Volume 



a 




exp 



p 



(27) 



where a is the lattice spacing. The upshot is that p, which is the fraction of plaquettes 
pierced by P-vortices (equals the probability that any given plaquette is pierced by a P- 
vortex) should scale like the string tension. 

Return now to the Creutz ratio plot in Fig. [7|, and in particular the data for the center- 
projected 1), where 



x(M) 



logW cp (l,l 



(2* 



where the "cp" subscript indicates that this is the center-projected Wilson loop. It is easy 
to see that 

W CP {1, l) = (l-p)+px (-1) = l-2p (29) 



so for small p (large (3) we have 



X(l,l)«2p 



(30) 



From the behavior of x(l, 1), which seems to (at least roughly) parallel the straight line 
shown, we see that p does appear to scale correctly. However, since fl3"0"| ) is approximate, 
it is better to plot the precise value of the P-vortex density 



P=^l-W cp (l,l)) 



(31) 



versus coupling (3, as shown in Fig. The straight line is the asymptotic freedom expression 
(last line of eq. ([2"7])), with the choice yj pj (6A 2 ) = 50. The scaling of P-vortex densities, 
at the larger (3 values, is rather compelling. There seems little doubt that P-vortices are 
locating physical, surface-like objects in the full Yang-Mills vacuum; objects which we have 
identified, in section 2.2 above, as center vortices. 
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Figure 9: Evidence for asymptotic scaling of the P- vortex density, defined as the fraction p 
of plaquettes pierced by P-vortices (one-sixth the average area occupied by P-vortices per 
unit lattice volume). The solid line is the asymptotic freedom prediction of eq. ([STD, with 

= 50. 



constant J pj (6A 



3 Cooled Vortices 



It has been argued persuasively by Teper |16] that the lattice cooling procedure can never, 
in a finite number of cooling steps, remove the asymptotic string tension extracted from 
sufficiently large Wilson loops. However, as the number of cooling steps increases, the 
area-law falloff sets in at increasingly large loop sizes; this means that for a lattice of any 
fixed volume, confinement is eventually lost. We would like to understand, in the vortex 
picture, how it is that the area law is lost for smaller loops while being preserved for larger 
loops, after a finite number of cooling steps. 

To answer this question, we need to know what happens to vortices as the lattice 
is cooled. As before, we use P-vortices found on the uncooled lattice to locate the center 
vortices, and to count the number of times n that vortices pierce a given lattice surface. The 
lattice is then cooled, using the constrained cooling procedure of Campostrini et al. JT7[ , 
and we can study what has happened to the configurations identified on the uncooled 
lattice. The first quantity of interest is Wi/Wq, where the Wilson loops are evaluated on 
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Figure 10: Variation of the ratio W\/Wq with the number of cooling steps. 



subensembles of the cooled, unprojected lattice, and the P-vortices are identified on the 
uncooled lattice. The result, from 0-20 cooling steps aX, (3 = 2.3, is shown in Fig. [TO] (all 
data in this section was obtained on a 16 4 lattice). A rough guide to the thickness of a 
vortex is the loop size for which Wi/Wq ~ 0. According to Fig. [TO], this happens for 2x3 
loops at cooling step 0, 3 x 4 loops at cooling step 5, and 4x5 loops at cooling step 10. 
The simplest interpretation is that the vortices become thicker as the cooling proceeds, 
with Wi/W reaching its asymptotic value at ever larger distance scales. 

The thickening of vortices with cooling explains how the area law is lost for smaller 
loops, but retained asymptotically for sufficiently large loops. The asymptotic string ten- 
sion is only obtained for loops whose dimensions are significantly larger than the vortex 
thickness; there is no area law falloff for loops whose size is very much smaller than the 
vortex thickness.^ Thus, as cooling begins, loops whose size is comparable to the vortex 
thickness lose their area-law falloff while the string tension of larger loops is unchanged. As 
cooling proceeds, vortex thickness increases, and the area law is lost for still bigger loops. 
However, after any number of cooling steps, there will always be loops (on a sufficiently 
large lattice) whose extension is large compared to the cooled vortex thickness, and whose 
asymptotic string tension is untouched. 

This picture is supported by the data for xo(R,R) vs - x{RiR) after 5 cooling steps 
(Fig. [11]), and after 10 cooling steps (Fig. [12]), again at (3 — 2.3. These figures should be 
compared with Fig. |5] above, which shows the same quantities on an uncooled lattice. We 
notice, particularly after 10 cooling steps, that the Creutz ratios for small loops have been 
drastically reduced. However, as R increases, the standard Creutz ratio x(-^> R) rises and 
seems to level out near the usual value of the asymptotic string tension at — 2.3. On the 

4 This point is discussed in much more detail in ref. M. 
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The Effect of P-Vortices at (3 = 2.3: Cool Step 5 



K 0.15 - 




Full Loops 
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Figure 11: Zero vortex Creutz ratio Xo(Rj R) an d the full Creutz ratio x(Ri R) vs - R> after 
5 cooling steps. 



other hand, the zero- vortex Creutz ratio Xo{R, R) is again tending to zero for large loops. 
At cooling steps (Fig. H), Xo(R, R) ~ at R — 5. At 5 and 10 cooling steps, Xo(5, 5) > 0, 
although the trend towards zero at increasing loop size is clear. A rough guess is that 
Xo(R, R)&0 for R such that Wx/Wq » -1. 

The message of figures p!0| - [T2| is that the vortices are still present on cooled lattices, 
and are still essential to confinement. However, the asymptotic values W\/Wq — > — 1 and 
Xo — > are obtained only at increasingly large loop area, as the number of cooling steps 
increases. This behavior, as well as the loss of area law falloff for smaller loops with cooling, 
seems to be nicely explained as being due to the "thickening" of the vortex core (which 
is the region of the center vortex which cannot be represented by a gauge transformation 
with discontinuity (|16|)). 

As just explained, our strategy is to locate the configurations of interest on the uncooled 
lattice, and then study what happens to these configurations as the lattice is cooled. But it 
is also interesting to ask whether our procedure for finding the center vortices, i.e. maximal 
center gauge combined with center projection, also works on the cooled lattices. The answer 
is "no." In Fig. |TB| we show how the P- vortex density p of eq. (^7j) falls drastically with 



cooling, if the P-vortices are identified by gauge-fixing and center-projecting the cooled 
lattice. A corresponding falloff is found in the Creutz ratios Xc P (R, R) of center-projected 
loops, with increasing cooling step, as seen, e.g., in Fig. |14] for R = 5. There is no trace of 
the typical "plateau" in string tension over some finite number of cooling steps, found in 
plotting the usual Creutz ratios. We have also computed Xc P (R, R) on lattices that have 
been "smoothed" according to the procedure in ref. |18| , and find that the center-projected 
Creutz ratios are reduced by about a factor of three on the smoothed lattices, as compared 
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The Effect of P-Vortices at ,3 = 2.3: Cool Step 10 
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Figure 12: Same as the previous figure, after 10 cooling steps. 



to the original lattices.^ Very similar phenomena have been reported in maximal abelian 
gauge, in refs. fl9j,p0"|. 



We stress, however, that the center vortices themselves, on a large lattice, are not 
removed by cooling. That fact seems evident in Fig. [1(HL2L What is lost on the cooled 



lattice is the efficacy of center-projection in finding all of the vortices. This suggests that 
our procedure for locating center vortices is sensitive not only to long-range fluctuations, 
but also to short-range features of some kind that are associated with these vortices. It 
would be interesting to know what these short-range features are.Q This question is under 
investigation, and we hope to return to it in a future publication. 



4 Gribov Copies 

The over-relaxation method of gauge- fixing, described in section 2, is not guaranteed to 
find the absolute maximum of the quantity R in eq. (0J); in general it will only find a local 
maximum. This is the well known "Gribov problem," which also afflicts the Coulomb, 
Landau, and maximal abelian gauges. It was in order to alleviate the problem that we 
have made three random gauge copies of each configuration used for data taking, and 
gauge-fixed each to obtain three "Gribov copies." We then used the Gribov copy with the 
largest value of R. 

Since Gribov copies of a given lattice configuration are not identical, it is interesting 
to study by how much, on average, they differ. In particular, to what extent are the 

5 We thank T. Kovacs for kindly supplying us with 100 smoothed lattice configurations. 
6 One clue is that the unprojected plaquette energy at the location of P-vortices, in the uncooled lattice, 
is significantly higher than the average plaquette energy. 
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p = [1-W (1,1)]/2 vs. Cool Step 
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Figure 13: The drop in P-vortex density p, identified by center-projection in maximal 
center gauge, on cooled lattices. 
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Figure 14: The drop in the center-projected Creutz ratio Xcp(5,5), corresponding to the 
drop in p, with cooling step. 
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positions of P-vortices correlated from copy to copy? If there is no correlation, then we 
can hardly rely on P-vortices to locate physical objects (i.e. the thick center vortices). On 
the other hand, some variation in the position of P-vortices, from copy to copy, should be 
permissible. Center vortices are rather thick, extended objects, and the precise "middle" 
surface of such configurations, which P-vortices are supposed to locate, may be somewhat 
ill-defined. 

To investigate quantitatively the correlation of vortices in different Gribov copies, we do 
the following: A Monte Carlo simulation is run at a given (3, taking data every 100 sweeps, 
and making 4 copies of each configuration chosen for data taking. The two "best" Gribov 
copies (in the sense of having the largest R values) are center-projected, and we denote the 
projected configurations by Z'(x) and Z"(x), with corresponding center-projected Wilson 
loops 

VV'(C) = Z'Z'...Z' 

W'{G) = Z"Z"...Z" (32) 
Then we compute the expectation value of the loop product 

< W'(C)W"(C) > (33) 
If the correlation of P-vortices in the two Gribov copies were perfect, then 

W{C) = W"(C) = ±1 (34) 

and therefore 

< W(C)W'(C) >= 1 perfect correlation (35) 

At the other extreme, if there were no correlation at all between the P-vortex positions in 
the two Gribov copies, then 

< W{C)W{C) > = < W{C) >< W'(C) > 

= exp[— 2crArea(C)] no correlation (36) 

where a is the string tension of the center-projected loops (same as the asymptotic tension 
of the unprojected loops). 

Small variations in P-vortex position among different Gribov copies are most likely to 
lead to a perimeter law falloff of the loop product (|33"D, at least in the limit of large loop 
area. Consider, for example, the following simple model: Take an / x J Wilson planar 
loop and assign to each of the plaquettes in its plane, both inside and outside the loop, 
a value +1 with probability (1 — /), and —1 with probability /. Each such configuration 
is supposed to represent a particular "Gribov copy" of center-projected plaquettes in the 
plane, and the value of the Wilson loop is WL(J, J) = (— l) n , where n' is the number of 
negative plaquettes in the minimal area. It is not hard to see, since the plaquettes are 
assumed to be uncorrelated, that, averaging over many configurations gives 

W' cp {I,J) = <W{I,J)> 

= (l-2/) /J = exp[-a/J] (37) 
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where a = — log(l — 2f) (the assumption that nearby plaquettes are completely uncor- 
rected is the main unrealistic feature of this model). From a given configuration of ±1 
plaquettes, we construct a second "Gribov copy" by allowing negative plaquettes to change 
their position by, at most, one lattice spacing. Then only changes in position of negative 
plaquettes bordering the loop perimeter can cause the loop product to differ from +1. 
Assign to a negative plaquette on the perimeter a probability q to cross into the loop, if 
it were outside, or to cross outside the loop, if it were inside. [| The new value of the loop 
is (— l) n , where n" is the number of negative plaquettes inside the loop in the second 
"Gribov copy." One then finds that in this model, defining A" = 4(1 + J), 

N y\n 

<w(i,j)w"(i,j)> = E(i-2g) n r(i-/) Ar - n - T7 T7— 

n=o n.{i\ n). 

= (1 - 2qfY {I+J) (38) 

which is a perimeter-law falloff. 

The above argument should also go through if the P-vortex positions vary by more 
than one lattice spacing among Gribov copies, so long as the variation is small compared 
to the size of the loop. If the variation in P-vortex position is comparable to the thickness 
of the center vortex, then our best chance to see perimeter-law falloff in the loop product 
(|33|), for comparatively small-size loops, will be at smaller values of coupling (3, where the 
vortex is relatively thin in lattice units. We have therefore chosen to do our simulation at 
a value of f3 = 2.1 which is just past the strong-to-weak coupling crossover. 



Figures [15] and [16| show our Monte Carlo results at (3 = 2.1 on a 14 4 lattice, with data 
from 400 configurations separated by 100 sweeps between configurations. Triangles show 
the data for loop products < W(C)W'(C) >, plotted vs. loop area, while the crosses 
are the values for no correlation, i.e. < W'(C) >< W"(C) >. The loop products are 
clearly far above the uncorrelated value and, from Fig. [15], do not seem consistent with an 
area-law falloff. In Fig. [IB] the loop product is plotted vs. loop perimeter. The straight 
line is drawn, somewhat arbitrarily, through data points at perimeter = 10, 18. It appears 
that the falloff in the loop product with perimeter is quite compatible with perimeter-law 
falloff, as predicted in our simple model. 

These results indicate that the variation in P-vortex position among different Gribov 
copies is relatively small - perhaps on the order of the vortex thickness, although we have 
not quantified this - and leads asymptotically only to a perimeter-law falloff for the loop 
product < W'(C)W"(C) >, indicating a strong correlation among Gribov copies. 



5 First Results in SU(3) 

All results presented in the previous sections support the idea that thick Z 2 vortices are the 
configurations dominating the SU(2) Yang-Mills vacuum. However, the vortex mechanism 

7 Plaquettes touching the corners of the loop should be treated a little differently from the other pla- 
quettes along the perimeter, but this is an inessential complication of the model, which we will ignore. 
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Correlation of Gribov Copies, Beta = 2.1 
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Figure 15: Expectation value of products of center-projected Wilson loops (triangles), 
evaluated in different Gribov copies, plotted vs. Loop Area. Crosses indicate the value for 
no correlation. 
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Figure 16: Expectation value of products of center-projected Wilson loops, evaluated in 
different Gribov copies, plotted vs. Loop Perimeter. 
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should not be specific to the SU(2) gauge group. In nature quarks appear in three colors, 
so a very urgent question is whether the observed phenomena survive the transition from 
SU(2) to 577(3). 

The maximal center gauge in SU(3) gauge theory is defined as the gauge which brings 
link variables U as close as possible to elements of its center Z3 = {e~ 2i7r / 3 /, 7, e 2 * 71 ' 3 !}. 
This can be achieved e.g. by maximizing the quantity 

i? = ^Re([Trf/ M (x)] 3 ), (39) 

X,fl 

or 

R! = £ I Tr U^x) | 2 . (40) 

X,fl 

We will concentrate here on the first choice. Center projection then amounts to replacing 
full link variable U^(x) by Z^(x), the closest center element. The residual unfixed local 
gauge symmetry is that of Z s . 

Fixing to the maximal center gauge in SU(3) gauge theory turns out to be much more 
difficult and computationally intensive than in the case of SU(2). The reason is that we 
have not succeeded in reducing the maximization to an underlying linear algebra problem as 
in SU(2) (see Section 2.1). We thus resorted to the method of simulated annealing |pl| , p2| , 
which was used for maximal abelian gauge fixing by Bali et al. [S3]. However, this method 



of maximal center gauge fixing converges to the maximum of R, Eq. (p9|), very slowly, 
which has forced us thus far to restrict simulations to small lattice sizes and to strong 
coupling. Tests of a more efficient maximization procedure are in progress. 

Before discussing the strong coupling results for SU(3), let us first show analogous data 
from SU (2) gauge theory. In Figure |17] we plot values of center-projected Wilson loops 
W(I, J) in maximal center gauge for (3 < 2.5. Broken lines connecting data points are 
just meant to guide the eye. Solid lines represent result of the lowest-order strong-coupling 
expansion (for unprojected loops) 

W(I,J)=(^j . (41) 

Monte Carlo data for projected loops agree with the lowest-order strong coupling expansion 
up to about (3 = 1.5. 

Our first results from 577(3) lattice gauge theory simulations come from an 8 4 lattice, 



for (3 values 1.0, 2.0, 3.0, 4.0, 5.0, 5.4 and 5.6. Figure |1S| shows center-projected Wilson 
loops together with the standard strong-coupling expansion to leading and next-to-leading 
order: 

IJ / T T _ \ 

(42) 



The data agree with lowest-order strong-coupling expansion up to (3 ~ 2; when next-to- 
leading term is taken into account, the agreement extends up to (3 — 4. 
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Figure 17: Center-projected Wilson loops vs. the strong-coupling expansion (solid lines) 
in SU (2) lattice gauge theory 




Figure 18: Center-projected Wilson loops vs. the strong-coupling expansion (solid lines) 
in SU (3) lattice gauge theory 
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Qualitatively, the situation at strong coupling looks much the same in SU (2) and 
SU(3): in both cases full Wilson loops are well reproduced by those constructed from 
center elements alone in maximal center gauge. Thus, center dominance is seen also in 
SU (3) gauge theory at strong coupling. 

An immediate task for the near future is to repeat our investigation of center dominance 
and the role of vortices in SU(3) lattice gauge theory for couplings in the scaling region. 
An absolutely crucial check of the validity of the vortex mechanism is that the evidence 
for vortices found in the SU(2) lattice theory is also found for the SU (3) gauge group. 

6 Summary 

It may be worth summarizing the results reported here: 

• P-Vortices locate center vortices. Vortex excitations in the center-projected 
configurations, in direct maximal center gauge, locate center vortices in the full, 
unprojected lattice. The evidence for this comes from the fact that 
W n (C)/W (C) -> (-l) n , and W odd (C) -> -W evn (C) in the limit of large loop area. 

• No vortices no confinement. When Wilson loops in SU(2) gauge theory are 
evaluated in subensembles of configurations with no vortices (or only an even number 
of vortices) piercing the loop, the string tension disappears. 

• Vortices, by themselves, account for the full string tension. The string 
tension of the vortex contribution to Wilson loops is found to match, quite accurately, 
the asymptotic string tension extracted from the full Wilson loops. 

• Vortex density scales. The variation of P- vortex density with coupling j3 goes as 
expected for a physical quantity with dimensions of inverse area. This is additional 
evidence that P-vortices locate physical, surface-like objects (center vortices) in the 
Yang-Mills vacuum (see also ref. ||15|| ). 

• Center vortices thicken as the lattice cools. This enables us to explain how 
the area law falloff is lost, after a finite number of cooling steps, for smaller loops, 
while the string tension remains unchanged for sufficiently large loops. 

• P-vortex locations are correlated among Gribov copies. There appears to 
be only modest sensitivity in P-vortex location to the choice of Gribov copy. 

• SU(3). There is preliminary evidence, on small lattices and strong couplings, of 
center dominance also in SU (3) lattice gauge theory. 

It is also worth mentioning some other results reported in refs. 0,0: 
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• Monopole loops lie on P-vortices p|. Monopoles, identified in the maximal 
abelian gauge, lie along center vortices, found in the indirect maximal center gauge, in 
a monopole-antimonopole chain. The non-abelian field strength of monopole cubes, 
above the lattice average, is directed almost entirely along the associated center 
vortices. Monopoles appear to be rather undistinguished regions of vortices, and 
may simply be artifacts of the abelian projection, as explained in ref. 0. 

• Center vortices are compatible with Casimir scaling [§. The "Casimir 
scaling" of the string tension of higher representation Wilson loops, at intermediate 
distance scales, has long been considered incompatible with the center vortex theory. 
Very recently, however, it has been argued that Casimir scaling is explained in terms 
of center vortices, if we take into account the fact that center vortices, unlike P- 
vortices, have a thickness which may be much greater than one lattice spacing.^ 

These results support the view that center vortices are responsible for quark confine- 
ment. 
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